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Ph! Abstract 



We consider the production of heavy quark jets at the Z-pole at the next-to- 
^ ■ leading order (NLO) using the Cambridge jet- algorithm. We study the effects 

of the quark mass in two- and three-jet observables and the uncertainty due 
to unknown higher order corrections as well as due to fragmentation. We 
found that the three-jet observable has remarkably small NLO corrections, 
which are stable with respect to the change of the renormalization scale, 
when expressed in terms of the running quark mass at the m^-scale. The size 
of the hadronization uncertainty for this observable remains reasonably small 
and is very stable with respect to changes in the jet resolution parameter y c . 
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I. INTRODUCTION 



During last few years a significant progress has been done in the understanding of the 
heavy quark jet-production in e + e~-annihilation both experimentally [] and theoretically. 

The DELPHI collaboration has measured the bottom-quark mass analyzing the 
e + e~-annihilation into the three-jet final state with heavy quarks using recent next-to- 
leading order theoretical predictions for this process [PHTOll- The DELPHI result for the 



Durham [IT jet clustering algorithm 

m b (m z ) = 2.67 ± 0.25(stat.) ± 0.34(W.) ± 0. 27 (theo.)GeV , (1) 

was the first measurement of the 6-quark mass far above the production threshold and it is 
the first experimental evidence (at the 2-3 sigma level) of the running of a fermion mass, 
as predicted by the Standard Model. Recently, the SLD collaboration has also analyzed its 
three- and four-jet data using Durham and several Jade-like jet algorithms JL2|. The value 
of the b-quark mass [|IJ| obtained from these data is compatible with the above DELPHI 
result. 

The heavy quark mass measurement was done under the assumption of the flavor in- 
dependence of the strong interactions with the value of the strong coupling constant, a s , 
fixed to its world average measured in other experiments. On the other hand, assuming a 
given 6-quark mass value obtained from low-energy measurements, and comparing the value 
of a s measured from the heavy quark three-jet final state with the one measured from the 
production of light quarks, one can perform a test of the flavor universality of the strong 
interaction. Such a test was performed recently 0,|l^,|l4[] and no deviation from the QCD 
prediction was found. The next-to-leading order QCD predictions with heavy quark mass 
corrections from references |§-10] were used in these studies. 



There are three main sources of uncertainties in the DELPHI analysis. The first one has 
statistical nature. The second error is due to the uncertainty in the hadronization correc- 
tions. It was evaluated || using different Monte-Carlo models simulating the hadronization 
process. The third one is due to our ignorance of higher order perturbative corrections in the 
theoretical predictions at the partonic level. The last uncertainty was estimated by varying 
the renormalization scale in the calculations and by using different renormalization schemes, 
i.e. expressing intermediate results in terms of either the perturbative pole quark mass or 
the running quark mass. 

The value of the 6-quark mass measured at the Z peak (|lj) is found to be in good 
agreement with the determinations of the fe-quark mass at low-energy from T- and 5-mesons 
spectroscopy ||15|| , when compared at the same scale. However, the uncertainties in m^mz) 
are larger. Thus, it would be desirable to reduce this error by finding new observables which 
may show a better theoretical and hadronization properties. 



1 See |TJ for a review of recent experimental results. 

2 The MS definition for the running mass at the m^-scale was used. 
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In this paper we study quark mass effects in heavy quark jet production by using the 
Cambridge jet clustering algorithm 
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We also study the possibility to reduce the 
uncertainties in the measurement of the 6-quark mass at the Z-pole. We consider two jet 
observables and estimate the errors in their theoretical predictions due to the unknown 
higher orders by varying the renormalization scale and considering different renormalization 
schemes. We discuss also the size of the uncertainty due to the hadronization process. Some 
preliminary results of this study were reported in [T3JT^ • 



II. THE CAMBRIDGE ALGORITHM AND THE DECAY Z -> 3JETS AT THE 

NEXT-TO-LEADING ORDER 



The Cambridge [0 jet clustering algorithm is a modified version of the popular 
Durham [[II]] algorithm that has been introduced recently in order to reduce at low y c the 
formation of spurious jets with low transverse momentum particles. Consequently, compared 
to Durham, it allows to explore regions of smaller y c while still keeping higher order cor- 
rections relatively small. It is important to note that at low y c the statistical experimental 
error for three-jet and four-jet production is expected to be smaller, and the sensitivity to 
the quark mass increases. 

In the Durham algorithm one finds the minimal test variable y^ defined as 

mini Eh E?) 

Vij = 2 ^-i^Z(l-cos%) , (2) 

for all possible pair-combinations of the particles and compare it with the jet-resolution 
parameter, y c . In eq. (0) Ei and Ej denote the energies of particles f\i and j, 9^ is the angle 
between their three-momenta and s is the center-of-mass energy squared. If y^ < y c , the 
two particles % and j are combined into a new pseudoparticle with momentum 

Pk=Pi+ Pj ■ (3) 

The procedure is repeated again and again until yij > y c for all pairs of (pseudo)particles. 
The number of (pseudo)particles at the end defines the number of jets. 

The Cambridge algorithm is defined by the same test variable, eq. (0), and the same 
recombination rule, eq. (|3|), as Durham. The new ingredient of the Cambridge algorithm is 
the so-called ordering variable 

Vij = 2(1 - cos%) . (4) 

In this algorithm, first, the pair of particles, which has minimal ordering variable v^, is 
selected. Then one computes y^ for this pair of particles, and, if y^ < y c , the two particles 
are recombined into a pseudoparticle according to eq. @. But if y^ > y c , the softer particle 
from this pair is assigned to a resolved jet. This last step is called "soft-freezing". 



3 By the word "particles" we mean here both the real hadrons detected at experiment and the 
partons entering the theoretical calculation. 
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Because of the additional step in the jet finding iterative procedure, the Cambridge 
scheme turns out to be more complex and has a number of peculiar properties [^0). Let 
us mention only one example. In Jade-like algorithms, including Durham, one can always 
define a transition value of y c , such that a multi-parton event classified as a n-jet event 
becomes a n + 1-jet event, when the value of y c is slightly decreased. However, as pointed 
out in [p0|| , this property is lost in the Cambridge algorithm, since, due to the presence of 
the ordering parameter, the sequence of clustering depends on the value of y c . As a result 
the number of jets is not a monotonic function of y c and it can change by a non unit number 
at some transition value of y c . 

With the above definitions one can show that the cross section of the e + e~-annihilation 
into three jets calculated at the leading order (LO) is the same in both Cambridge and 
Durham algorithms. This happens because only three-parton final state configurations con- 
tribute at LO to the three-jet cross section. Instead, the four-jet production cross section at 
the LO is different in the two schemes. 

At the next-to-leading order (NLO) the predictions for the three-jet production cross sec- 
tion for the two algorithms are different. Schematically, the NLO calculation of e + e _ —> 3jets 
was performed as follows. In this case the three-jet cross section receives contributions not 
only from one-loop corrected three-parton final states, but also from four-parton processes. 
In the latter process two of the four partons are combined in order to produce a three-jet final 
state. The ultraviolet (UV) divergences encountered in the calculation of the three-parton 
contribution at the one-loop level were removed by the renormalization of the parameters 
of the QCD Lagrangian. The infrared (IR) divergences [] remaining in this part, which are 
due to the presence of massless gluons in the loop, were canceled in the final result for the 
three-jet transition probability by adding an appropriate contribution from the four-parton 
final state. In the latter contribution, which is a purely tree-level one, the IR divergences 
appear due to the radiation of soft or/and collinear massless gluons. To separate the IR 
divergent part of the four-parton contribution, the phase-space slicing method (see |H] and 



references cited therein) has been used. In this method the integration over a thin slice at 
the edge of the phase-space (containing the soft and the collinear singularities) is performed 
analytically. Then, the IR singularities coming from three- and four-parton final states are 
canceled analytically. The remaining finite pieces from both three-parton and four-parton 
processes are integrated numerically over the three-jet phase-space defined by the specific jet 
algorithm. The four-jet cross section at the leading order, which is IR finite, is also obtained 
by numerical integration over the four-jet part of the four-parton phase-space. 

Details of the NLO calculation for the Durham and some other popular jet clustering 
algorithms were presented in In the case of the Cambridge algorithm, although all 

principal calculational steps of the three-jet heavy quark production in e + e~-annihilation at 
NLO remain the same as in other algorithms, the practical implementation of this scheme 
turned out to be more involved due to more complex realization of the Cambridge jet finder. 



4 Dimensional regularization was used to regularize both ultraviolet (UV) and infrared (IR) sin- 
gularities in the whole calculation. 



4 



III. THE OBSERVABLES 



In this paper we study in detail the following ratio of three-jet rates in the Cambridge 
jet-algorithm 

r>b£ _ r 3j(^)/r b 

* 3 - r^/i* • (5) 

In the above equation r| ■ and T b stand, respectively, for the three-jet and the total decay 
widths of the Z-boson with a 6-quark in the final state. Analogously, the quantities with 
the superscript I denote the sum of the decay widths into light quarks (i = u, d, s) which all 
are considered massless. 

The ratio Kg can be written in the form of the following expansion in a s 

Rf = 1 + ^a (y c ) + r b (b (r b ,y c ) + ?M h^yS] , (6) 

where r b = M^/s, with M b the heavy quark pole mass, and s = m? z at the Z peak. 

Let us remark that the double ratio in eq. (|5|) differs slightly from the one, R^ d , considered 
22|,|4| with a normalization to the Z-decay width of only one light flavor, the <i-quark. 



m 



In contrast to R z , such double ratio is equal to unit for a vanishing 6-quark mass. The 
main difference between the two observables, _R[f and R^, is due to the triangle one-loop 
diagrams II, which give a non-zero contribution even in the case of massless 6-quarks 



taken into account by a function ao(y c ) in eq. ([]). The difference is, however, very small 
numerically, smaller than a 0.2%. 

The functions &o and b± in eq. (Q) describe, accordingly, the quark mass effects at the 
leading and the next-to- leading order in the strong coupling and depend on the jet clustering 
scheme. Although, for convenience, the leading polynomial dependence on has been 
factorized out in eq. @, the exact dependence on the heavy quark mass is kept in the 

functions b (r b ,y c ) and &i(r 6 ,y c ). 

Using the known relationship between the perturbative pole mass and the MS scheme 
running mass 
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M l h = m 2 M 



l 2a e (ji) 




(7) 



we can re-express eq. @ in terms of the running mass m(,(/x). Then, keeping only terms of 
order 0(a s ) we obtain 



R 3 =l + — a Q (y c ) + r b (fi) \b {r b ,y c ) + — h{r b , y c , fi) ) , (8) 

7T y 7T J 

where f b {p) = ml(/J,) / m| and the new function b\ is related to bo and &i, introduced in eq. (H, 
via 



bi(r b ,y c ,fi) = h(r b ,y c ) + 2b (r b ,y c ) ( ^ - logr 6 + log ) • ( 9 ) 

6 Tflr. 
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Effectively, the use of m&(/x), instead of Mb, corresponds to the use of a different renormal- 
ization scheme. Although at the perturbative level both expressions, eq. (|]) and eq. 
are equivalent, numerically they give different answers since different higher order contribu- 
tions are neglected. The spread of the results gives an estimate of the size of higher order 
corrections. 

As discussed in the previous section, the phase-space integration (up to five-fold) in the 
calculation of the NLO decay width of the Z-boson into three jets is done numerically. This 
numerical integration is rather time consuming. Hence, we found it very convenient to fit the 
numerical results with relatively simple analytical functions. Because very small y c -values 
are considered in the case of the Cambridge scheme, these fits are more complex and involve 
more parameters than the ones for the Durham algorithm described in Q. A Fortran code 
containing the fits to the functions 60 and b\ (or b\ )f] can be obtained from the authors upon 
request. The numerical results for R$ are presented in the next section. 

In the next section we also give numerical results for the ratio of differential two-jet rates, 
defined as follows 

D u _ \ V UVc + Ay c /2) ~ T%(y c - Ay c /2)}/T b 
2 [Y%(y c + Ay e /2) - V%{y c - ' 1 } 

Here, r| • and T^j denote the two-jet decay widths of the Z-boson with a 6-quark and light 
quarks in the final state, correspondingly. The two-jet decay width at the order 0(a 2 s ) is 
calculated from the three- and the four-jet widths through the identity 

where q is the quark flavor and r|- is the four-jet decay width at the leading order. The 
value of Ay c in eq. (0) should be chosen small enough. We fix Ay c = 0.001 in the numerical 
analysis. 

The differential ratio D?f is interesting because it contains different information than the 
ratio i?! . In addition, while values of B% measured at different y c are strongly correlated, 
the differential rate U% can be analyzed as a function of y c . The whole consideration of 
i?3 discussed above is also applied here. For D% we use expansions in a s similar to those 
in eq. @ and eq. (||), see 0, and we fit the corresponding LO and NLO numerical results 
to simple analytical functions equivalent to bo and b\ (t>i). 

IV. PERTURBATIVE RESULTS 

In fig. m we present the results for the two observables studied, and Dtf, as functions 
of the jet-resolution parameter y c in the Cambridge algorithm. We plot the NLO results 
written either in terms of the pole mass (eq. (|Bp) with Mb = 4.6 GeV, or in terms of the 
running quark mass at rriz (eq. (§)) with rribimz) = 2.8 GeV. The renormalization scale is 



5 Although, the two functions b\ and b\ are related via eq. @ we performed independent fits for 
these functions. 
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fixed to ji = rriz and a s (rriz) = 0.118. For comparison we also show -Rjf and Dtf at LO 
when the value of the pole mass, Mf,, or the running mass at mz, m^mz), is used for the 
quark mass. Note that one can not distinguish between different definitions of the quark 
mass in the lower order calculation. Mass effects monotonically grow for decreasing y c , they 
are very significant for both observables and in the case of exceed 10% for small values 
of y c . 

From this figure one sees a remarkable feature of the NLO result in the considered range 
of y c , 0.005 < y c < 0.025, for the Cambridge scheme : the NLO corrections, in the case when 
the running mass is used are significantly smaller, especially for R$ , than the corrections 
in the case with the parameterization in terms of the pole mass. In other words, using the 
running mass at the m^-scale in the LO calculations takes into account the bulk of the NLO 
corrections. This situation, although does not guarantee, suggests that also next-to-next-to- 
leading and higher order corrections are small for the observables parameterized in terms of 
the running mass at the m^-scale, i.e. one has a better description of mass effects in terms 
of a short distance parameter, m^mz), than in terms of a low-energy parameter like the 
perturbative pole mass. 

The theoretical prediction for the observables studied contains a residual dependence on 
the renormalization scale fi: when written in terms of the pole mass it only comes from the 
/^-dependence in a s (/i), when written in terms of the running mass it comes from both a s (/t) 
and the incomplete cancellation of the /^-dependences between mf,(/i) and the logs of /i which 
appear in eq. The dependence on /i is usually regarded as an estimate of the effect of 
the unknown higher order perturbative corrections. In fig. |2|.a we present the /t-dependence 
of the two NLO predictions, the pole mass prediction (NLO-Mb) given by eq. (^|) and the 
running mass prediction (NLO-m^m^)) given by eq. (||]), for the ratio i?§ in the range 
m z /l0 < n < m z at a fixed value of y c . We use the following one-loop evolution equations 

= a ^A , mM = mh[mz) k-io/Po j 



where a(/i) = a s (fi)/ii, K = 1 + a (mz)j3o log (/i 2 /m|) and 

7o = 1 



11 - -N F 
3 



with Np = 5 the number of active flavors, to obtain a s (fi) and mj(/i) from a s (mz) = 
0.118 and rrib(mz) = 2.8 GeV. The NLO-mb(mz) result (dashed line) shows a remarkable 
stability with respect to the variation of the renormalization scale and the corrections with 
respect to the LO prediction (LO-m^mz)) remain small for all the values of \i. Instead, the 
NLO-Mb prediction (dotted line) has noticeably stronger dependence on the renormalization 
scale. The NLO corrections in this case remain sizable for all the values of fi and increase 
for decreasing fi. Note also that, as one would expect, for low values of /i the two NLO 
predictions, in terms of the running mass and in terms of the pole mass, become very close 
to each other. 

For a given value of B!§ we can solve eq. (||) (or eq. (§)) with respect to the quark mass. 
The result, shown in fig. |.b for a fixed value of , depends on which equation was used 
and has a residual dependence on the renormalization scale /i. The curves in fig. 0.b are 
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obtained in the following way: first from eq. @ we directly obtain for an arbitrary value of 
\i between m z and m z /10 a value for the bottom-quark running mass at that scale, m&(/x), 
and then using eq. ([□]) we get a value for it at the Z-scale, rrib{m z ). Second, using eq. (||) 
we extract, also for an arbitrary value of /i between m z and m z /10, a value for the pole 
mass, Mb. Then we use eq. (0) at fi = Mb and again eq. (|11|) to perform the evolution from 
/i = Mb to fi = m z and finally get a value for m,b(m z ). The two procedures, we denote them 
NLO-mb{m z ) and NLO-M^ respectively, give different answers, since different higher orders 
have been neglected in the intermediate steps. The maximum spread of the two results in 
the whole /i-range under consideration can be interpreted as an estimate of the size of higher 
order corrections, i.e. of the theoretical error in the determination of the bottom-quark mass 
from the experimental measurement of i?f . 

We see from fig. ^.b that the first approach is very stable with respect to the choice of the 
scale used in eq. (H). The obtained 6-quark mass, m b (m z ), varies only ±50 MeV when the 
scale is varied in the range // = m z and \i = m z /10. In the same range of //, the estimated 
error in the Durham algorithm was found [[| to be ±200 MeV. On the contrary, if one uses 
eq. (^|) the extracted quark mass has a strong scale dependence, specially for small /i-values 
and the estimated error is very sensitive to the choice of the smallest possible value of the 
renormalization scale. Cutting as before at \i = m z /10, the extracted pole mass varies in 
the range ±300 MeV which is translated into ±240 MeV for mb(m z ). Let us note that a 
further ±20 MeV should be added due to the uncertainty in the strong coupling constant 
Aa s (m z ) = ±0.003. 

Although, our observables are formally of order 0(a s ) and, therefore, compatible with the 
use of one-loop renormalization group equations (RGE) to connect the running parameters 
at different scales, as a check of the stability of our results we have also repeated the analysis 
using two-loop evolution equations [^5J 

/ v a(m z ) 

K + a{m z )b x [h + a{m z )b x X ~ \ + L ) ' 

m b (m z ) = m b (n) K°° VV*^ ~ • ( 12 ) 

1 ± a(fjL) Ci 

where L = log if and c\ — g\ — big with b\ = Pi/fio, 9% = li/Po and 

ft = 1 U - ™N F ] , 7l = 1 m - ^AVI . (13) 
Hl 16 [ 3 J ' ' 16 L 3 9 J y ! 

The use of two-loop RGE corresponds to the dashed lines in Fig ||.b. Again a value of 
the quark mass extracted via the running mass parameterization remains more stable with 
respect to variation of the scale fi and changes only slightly. The mass extracted through 
the pole mass receives a significant shift of 200 MeV when the two-loop RGE are used. 

V. HADRONIZATION CORRECTIONS 

In the DELPHI analysis on the measurement of the tVquark mass effects based on the 
Durham jet clustering algorithm the impact of the fragmentation process on the observable 
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R*§ was studied and quantified by adding in quadrature two different source of errors. The 
first uncertainty, <J tU n, was obtained by varying the most relevant parameters of the string 
fragmentation model incorporated in JETSET J26] within an interval of ±2cr from its central 



value as tuned by DELPHI |27j and explained in reference [[| . The second uncertainty, <J mo d, 
was the result of analyzing the dependence on the fragmentation model itself by comparing 
the HERWIG [p^j model with JETSET . The difference on the fragmentation correction 
factors obtained for each model were considered as a source of systematic errors. This in 
fact was the largest contribution to the total error of the measurement. The final correction 
adopted was the average of those two models and the fragmentation model uncertainty 
{v-mod) w as taken to be half of their difference. The total error due to the lack of knowledge 
on the hadronization process was expressed as: 



OhadiVc) = \/ VtuniVc) + <jl od (yc) (14) 

which at y c = 0.02 in the Durham scheme was Uhadiyc) = 0.007 and its dependence as a 
function of y c is shown in Fig. ^. The decision of the measurable y c interval region, y c > 0.015 
was also connected to the fragmentation correction which was required to be relatively flat 
and the four-jet contribution small (< 2%). For comparison purposes, an equivalent analysis 
has been performed using the new Cambridge jet reconstruction algorithm and the results 
obtained are also presented in Fig. || A larger flat ?/ c -region is observed in the case of 
Cambridge with respect to Durham which can be extended up to y c = 0.004 with the four 
jet contribution still being small, <8%. The total absolute error is higher for Cambridge 
than for Durham but the relative sensitivity to the mass correction is higher for Cambridge 
at y c values around 0.005 than for Durham at y c values around 0.02. In Fig. [3] the difference 
between the theoretical prediction of Rf at LO in terms of the pole mass, = 4.6 GeV 
with respect to that obtained using the running mass m^imz) = 2.8 GeV is also shown. A 
higher sensitivity to this difference is again found for the new Cambridge jet algorithm in 
the valid, flat, ?/ c -region (y c > 0.004) thus enabling to test more significantly which of the 
two predictions agrees better with data. 



VI. CONCLUSIONS 

We have calculated the next-to-leading order QCD corrections to the heavy quark three- 
jet production cross section in e + e~-annihilation as well as the leading order four-jet pro- 
duction cross section using the new Cambridge jet clustering algorithm. The hadronization 
corrections were also estimated. Comparing with previous studies, this algorithm allows to 
extend the analysis into a region of smaller values of the jet resolution parameter, down to 
y c ~ 0.004, where the sensitivity to the heavy quark mass effects increase. 

In particular we have studied in detail the double ratio R h ^ and the differential double 
ratio D%- We have compared the NLO results expressed in terms of the perturbative pole 
mass and in terms of the running mass of the heavy quark at the mz scale. We found that 
the NLO corrections in the case when the running mass was used are remarkably small. This 
is especially true for Kg , where tree-level expressions in terms of mb(mz) give a very good 
approximation to the complete NLO result, which, when expressed in terms of m^imz)-, is 
almost independent of the renormalization scale. In contrast, the calculations done in terms 



9 



of the pole quark mass have sizable NLO corrections with a strong renormalization scale 
dependence. The hadronization corrections also favor the use of the Cambridge algorithm 
with respect to Durham by keeping a relatively stable uncertainty for small y c . 

Summarizing, the results of this paper indicate that a new determination of the 6-quark 
with the Cambridge jet algorithm will improve our present understanding on quark mass 
effects at LEP/SLC energies. 
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FIGURES 




FIG. 1. The observables K§ and as a function of y c in the Cambridge algorithm at the NLO. 
The dotted lines give the NLO corrected values using eq. (^) for a pole mass of M& = 4.6 GeV. The 
dashed lines give the observables at the NLO using eq. (||) for a running mass of mb(mz) = 2.8 GeV. 
The renormalization scale is fixed to \i = mz and a s {mz) = 0.118. For comparison we also plot 
the LO results for Mb = 4.6 GeV (lower solid lines) and rrib(mz) = 2.8 GeV (upper solid lines). 
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FIG. 2. a) Renormalization scale dependence for a fixed value of y c . Same labels as in Fig. |]. 
b) Extracted value of m^mz) from a fixed value of using either the pole mass expression 
(NLO-Mb) in eq. (0) or the running mass expression (NLO-m^mx)) in eq. (^) as explained in the 
text. Solid lines obtained by using one- loop running evolution equations to connect the results at 
different scales and dashed lines obtained by using two-loop expressions. 
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FIG. 3. Comparison of the hadronization uncertainty {a had) obtained when using either the 
Cambridge or Durham algorithm. The dashed curve shows the mass correction at LO for the 
pole mass, = 4.6 GeV. The dotted curve indicates the value of the difference between the 
LO predictions for the two mass values, M& = 4.6 GeV and mb(mz) = 2.8 GeV. The Cambridge 
algorithm is observed to have a larger stable region on y c than Durham reaching at the same time a 
higher sensitivity to both, the mass correction and the difference between the two LO predictions. 
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